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ABSTRACT 

The majority of galaxy mergers are expected to be minor mergers. The observational 
signatures of minor mergers are not well understood, thus there exist few constraints 
on the minor merger rate. This paper seeks to address this gap in our understanding 
by determining if and when minor mergers exhibit disturbed morphologies and how 
they differ from the morphology of major mergers. We simulate a series of unequal- 
mass moderate gas-fraction disc galaxy mergers. With the resulting g-band images, we 
determine how the time-scale for identifying galaxy mergers via projected separation 
and quantitative morphology (the Gini coefficient G, asymmetry A, and the second- 
order moment of the brightest 20% of the light M20) depends on the merger mass ratio, 
relative orientations and orbital parameters. We find that G~M2o is as sensitive to 9:1 
baryonic mass ratio mergers as 1:1 mergers, with observability time-scales ~ 0.2—0.4 
Gyr. In contrast, asymmetry finds mergers with baryonic mass ratios between 4:1 and 
1:1 (assuming local disc galaxy gas- fractions). Asymmetry time-scales for moderate 
gas- fraction major disc mergers are ~ 0.2—0.4 Gyr, and less than 0.06 Gyr for moderate 
gas- fraction minor mergers. The relative orientations and orbits have little effect on the 
time-scales for morphological disturbances. Observational studies of close pairs often 
select major mergers by choosing paired galaxies with similar luminosities and/or 
stellar masses. Therefore, the various ways of finding galaxy mergers (G — M20, A, 
close pairs) are sensitive to galaxy mergers of different mass ratios. By comparing 
the frequency of mergers selected by different techniques, one may place empirical 
constraints on the major and minor galaxy merger rates. 

Key words: galaxies:interactions - galaxies: structure - galaxies:evolution 
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1 INTRODUCTION 

The vast majority of galaxy mergers are expected to be 
mergers between unequal-mass systems. Within the stan- 
dard cosmological framework, the frequency of dark matter 
halo mergers increases with the mass ratio of the merger 
(e.g. Lacey & Cole 1993; Fakhouri & Ma 2008; Genel et 
al. 2009). 'Minor' mergers of galaxies with total mass ratios 
greater than 3:1 — 4:1 are more common than 'major' merg- 
ers with mass ratios closer to 1:1. Cosmological-scale simu- 
lations and semi-analytic models suggest that minor merg- 
ers play a significant role in triggering star-formation in the 
early universe (e.g. Somerville, Primack, & Faber 2001), ac- 
creting mass onto galaxies (e.g. Guo & White 2008; Genel et 
al. 2009; Stewart et al. 2009), and growing the stellar haloes 



* NOAO Leo Goldberg 
t W.M. Keck Fellow 



Fellow 



of massive galaxies (e.g. Johnston, Hernquist, & Bolte 1996; 
Bell et al. 2008). 

However, there exist few observational constraints on 
the minor merger rate and its evolution with time. Galaxy 
mergers may be identified either by searching for pre-merger 
systems of galaxies close in projected separation and relative 
velocity, or for galaxies with disturbed morphologies. Minor 
merger studies with close pairs are hampered by the high 
contrast ratio between the primary and satellite galaxies. It 
is challenging to obtain spectroscopic velocities or reliable 
photometric redshifts for the satellite galaxies with masses 
and luminosities ~ a tenth of the primary galaxies. Addi- 
tionally, the contamination rate from background galaxies 
increases dramatically as the luminosity difference between 
the primary and satellite galaxies increases. 

It has generally been assumed that the morphological 
effects of minor mergers are subtle. High resolution N-body 
simulations of gas-rich minor mergers find that such events 
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Figure 1. u — r — z composite colour images including dust extinction for the 3:1 baryonic mass ratio tilted prograde-prograde G3G2Pf 
simulation. The G3 primary galaxy is viewed face-on (camera 0). The images show the initial galaxies [1], the first pass [2], the maximal 
separation after the first pass [3], the second pass [4], the final merger [5], the post-merger ^ 0.5 Gyr after the merger [6], and the 
remnant ~ 1 Gyr after the merger [7]. The field of view for panels 1 and 3 is 200 kpc, while the field of view for the other images is 100 
kpc. 



may thicken but not destroy the discs of the primary galaxies 
(Moster et al. 2009; Purcell, Kanzantzidis, & Bullock 2009; 
Robertson et al. 2006; Abadi et al. 2003), as may minor 
mergers on retrograde or circular orbits (e.g. Velazquesz & 
White 1999; Hopkins et al. 2009a). Similarly, minor mergers 
may deposit material onto the outer regions of high redshift 
compact spheroids, but do not significantly change their 
inner density structure (Naab et al. 2009; Hopkins et al. 
2009b). A study of a 5:1 pure N-body galaxy merger sim- 
ulation found that the merger never showed asymmetries 
strong enough to meet the standard merger criterion (Con- 
selice 2006). However, recent studies of the visual classifi- 
cation of galaxy mergers suggests that ^ 50% of visually 
disturbed galaxies may result from minor mergers (Jogee et 
al 2009; Darg et al. 2009). Minor mergers and interactions 
may also produce lopsided discs (e.g. Rudnick & Rix 1998; 
Bournaud et al. 2005; Reichard et al. 2008). Therefore the 
morphological effects of minor mergers during the merger 
event may be more significant than previously thought. 

With a better physical understanding of why and when 
merging galaxies appear morphologically disturbed, we can 
constrain the observed galaxy minor and major merger rate. 
Most previous theoretical studies of minor mergers have fo- 
cused on the structures of the merger remnants (e.g. Hopkins 
et al. 2009; Naab et al. 2009; Robertson et al 2006; Bour- 
naud, Jog, & Combes et al. 2004; Johnston et al. 1996) or 
spanned a very limited range of parameter space (Conselice 
2006). Bournaud et al. (2005) performed a detailed study of 
the lopsidedness of disc galaxy mergers, and found that gas- 
rich minor mergers could appear strongly lopsided. In this 
work, we explore the strong morphological signatures pre- 



dicted by simulations of disc mergers with a range of mass 
ratios and orbital parameters. 

Translating the number of observed merger candidates 
into a merger rate requires the assumption of an observ- 
ability time-scale - the time during which one would have 
identified the system as merging. Until now, the time- 
scales for detected disturbed morphologies have been poorly 
constrained. The observability time-scale for a particular 
merger may depend on (1) the method used to identify the 
merger; (2) the merger parameters (mass ratio, gas prop- 
erties, bulge/disc ratio, orbits, dust content); and (3) the 
observational selection (observed wavelength, viewing an- 
gle, spatial resolution). Cosmological-scale numerical simu- 
lations currently do not have the spatial resolution to di- 
rectly determine the cosmologically-averaged observability 
time-scale for each method. Therefore, one is required to 
use a suite of galaxy-scale numerical simulations which span 
a large range of input merger parameters to constrain the 
observability time-scales for the different input parameters. 
Given a sufficiently broad range of merger parameters, the 
observability time-scales for each parameter set may then 
be weighted by the probability distribution of the mass ra- 
tios, gas fraction, etc. which can be computed from current 
cosmological-scale simulations. An additional complication 
is that galaxy-scale numerical simulations typically track the 
distribution of 'particles' (star, gas, and dark matter), as 
opposed to the projected fight distribution at a particular 
wavelength (which is what is observed). 

We continue the work first presented in Lotz et al. 
2008b, hereafter Paper 1. In Paper 1, we used GADGET 
N-body/SPH simulations processed through the radiative 
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transfer code SUNRISE to obtain the rest-frame g mor- 
phologies of equal-mass gas-rich disc mergers. We found that 
the time-scales for identifying a particular galaxy merger 
could be quite different depending on the gas-fraction of the 
merging galaxies and the method used to find the merger. 
For a given method, the time-scales associated with strong 
morphological disturbances are not a strong function of or- 
bit, orientation, total mass, or supernovae feedback. The ef- 
fects of numerical resolution, dust, viewing angle, and spatial 
resolution are also presented in Paper 1. 

In this paper, we focus on the effect of mass ratio on 
the quantitative morphologies and projected separations of 
simulated mergers of disc galaxies with gas fractions, bulge- 
to-disc ratios, and masses tuned to match local galaxies. 
In §2, we describe the simulations, the properties of the 
initial galaxies, and the range of merger parameters. In 
§3, we briefly describe the analysis of the resulting simu- 
lated g-band images and the criteria for identification as 
a merger by morphology and projected separation. In §4, 
we discuss the resulting observability time-scales and their 
dependence on the merger mass ratios, orbits, and relative 
orientation. We will present the effects of gas fraction in 
a companion paper (Paper 3; Lotz et al. 2009). Those fa- 
miliar with our approach from Paper 1 may skip to §4—5 
for the results and discussion. The simulated g band im- 
ages and morphology tables will be available in 2010 at 
the Muhimission Archive at STScI (MAST) as a High- 
Level Science Products (HLSP) contribution "Dusty In- 
teracting Galaxy GADGET-SUNRISE Simulations" (DIG- 
GSS); http://archive.stsci.edu/prepds/diggss . 



2 SIMULATIONS 

Here we briefly describe the galaxy merger simulations and 
initial conditions. 



2.1 GADGET N-Body/SPH simulations 

The details of these simulations, their global star-formation 
histories, and their remnant properties are discussed in Cox 
et. al (2004, 2006, 2008). AU of the simulations presented 
in this work were performed using the N-Body/SPH code 
GADGET (Springel, Yoshida, & White 2001). While we 
use the first version of GADGET (Springel et al. 2001), the 
smoothed particle hydrodynamics are upgraded to use the 
'conservative entropy' version that is described in Springel & 
Hernquist (2002). Each galaxy is initially modeled as a disc 
of stars and gas, a stellar bulge, and a dark matter halo, 
with the number of particles and masses for each compo- 
nent given in Table 1. The stellar and dark matter particles 
are coUisionless and are subject to only gravitational forces. 
The gas particles are also subject to hydro-dynamical forces. 
The baryonic and dark matter particles have gravitational 
softening lengths of 100 pc and 400 pc respectively. The SPH 
smoothing length for the gas particles indicates the size of 
the region over which the particle's hydrodynamic quanti- 
ties are averaged and is required to be greater than half 
the gravitational softening length or > 50 pc. The radiative 
cooling rate A„et (p, u) is computed for a primordial plasma 
as described in Katz et al. (1996). 



Gas particles are transformed into coUisionless star par- 
ticles assuming the Kennicutt-Schmidt law (Kennicutt 1998) 
where the star-formation rate depends on the local gas den- 
sity pgas ■ These new star particles have typical masses ~ 10^ 
M0, and are assigned ages based on their formation time 
and metallicities based on the metallicity of the gas particle 
from which they are spawned. We adopt the instantaneous 
recycling approximation for metal production whereby mas- 
sive stars are assumed to instantly become supernovae, and 
the metals produced are put back into the gas phase of the 
particle. In this version of GADGET, metals do not mix 
and remain in the gas particle in which they are formed. 
The enriched gas contribution from stellar winds and Type 
la supernovae are ignored. Unlike the metals, there is no 
recycling of hydrogen and helium to the gas. 

Feedback from supernovae is required to produce stable 
star-forming discs. We adopt a model in which the supernova 
feedback energy is dissipated on an 8 Myr time-scale, and 
have a equation of state P ~ p^aa- No active galactic nuclei 
(AGN) are included in these simulations. As we discussed in 
Paper 1, the exclusion of AGN feedback will not affect the 
morphological disturbance time-scales calculated here but 
may affect the appearance of the merger remnants. 

2.2 SUNRISE Monte Carlo radiative transfer 
processing 

SUNRISE is a parallel code which performs full Monte 
Carlo radiative-transfer calculations using an adaptive-mesh 
refinement grid (Jonsson 2006; Jonsson et al. 2006; Jonsson, 
Groves, & Cox 2009). We use SUNRISE v2 to create g-band 
images at least 30 timesteps for each merger simulation. For 
each GADGET simulation timestep, SUNRISE assigns a 
spectral energy distribution to each star particle using the 
STARBURST99 population synthesis models (Leitherer et 
al. 1999). 

The metallicities of the gas and stars of the initial 
galaxy models decline exponentially with the radius of the 
disc. The density of dust is linearly proportional to the den- 
sity of metals in the gas. The central metallicities and gra- 
dients scale with the mass of the galaxy. The details of the 
inital galaxy metallicities, gradients, and dust extinctions 
are given in Rocha et al. (2008). 

Given a particular simulation geometry and view- 
ing angle, SUNRISE v2 performs the Monte-Carlo radia- 
tive transfer calculation for 20 wavelengths from the far- 
ultraviolet to the mid-infrared and interpolates a resulting 
spectral energy distribution of 510 wavelengths including the 
effects of absorption and scattering. Images are created for 
eleven isotropically positioned viewpoints ("cameras"). In 
Figures [T] [2] and|31 we show examples of composite u~r ~ z 
images for the G3G2, G3G1, and G3G0 tilted prograde- 
prograde simulation viewed face-on (camera 0) at multiple 
timesteps. 

2.3 Initial Galaxy Models 

In order to sample the parameter space spanned by many 
present-day galaxies, we explored mergers between model 
galaxies with masses, bulge-to-disc ratios, and gas fractions 
motivated by SDSS estimates of typical local galaxies (Table 
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Figure 2. u — r — z composite colour images including dust extinction for the 9:1 baryonic mass ratio tilted prograde-prograde GSGlPt 
simulation. The viewing angles, merger stages and image scales are the same as Figure [T] 



Table 1. Initial Galaxy Conditions 



Model 


AT ^ 




M ■ ^ 

(M@) 




(Me) 


{MqT 


M*f 
(M~ 


(Mq) 


f h 




(kpc) 


(kpc) 


(kpc) 


G3 


2.4, [1.2,0.5,0.5,0.2] 


10^ 


1.2 • 10^2 


6 


6.2 ■ lO^" 


4.1 ■ 10^" 


8.9- lO'' 


1.2 ■ lO^O 


0.14 


0.19 


2.85 


0.62 


8.55 


G2 


1.5, [0.8,0.3,0.3,0.1] 


10^ 


5.1 • 10" 


9 


2.0 ■ lO^" 


1.4- 10^" 


1.5 ■ 10^ 


4.8 ■ 10^ 


0.08 


0.24 


1.91 


0.43 


5.73 


Gl 


9.5, [5.0,2.0,2.0,0.5] 


10-* 


2.0- 10" 


12 


7.0 ■ 10^ 


4.7 ■ 10^ 


3.0- 10** 


2.0- 10^ 


0.04 


0.29 


1.48 


0.33 


4.44 


GO 


5.1, [3.0, 1.0, 1.0,0.1] 


10-* 


5.1 ■ 10"' 


14 


1.6 ■ 10^ 


9.8 ■ 10* 


2.0- 10^ 


6.0- 10* 


0.01 


0.38 


1.12 


0.25 


3.36 



" Number of particles in total, dark matter, gas, stellar disc, and stellar bulge for GADGET simulations 
^ Virial mass 

Dark matter halo concentration 

Baryonic mass 
^ Mass of stellar disc 
f Mass of stellar bulge 
9 Mass of gaseous disc 
^ Fraction of baryons in the bulge 
* Fraction of baryons in gas 

Scalelength of stellar disc 
^ Scalelength of bulge 
' Scalelength of gaseous disc 



[T]). Each galaxy model contains a rotationally supported disc 
of gas and stars, a non-rotating stellar bulge, and a dark 
matter halo (Table [1} • A detailed description of the galaxy 
disc models can be found in Cox et al. (2006, 2008), Jonsson 
et al. (2006) and Rocha et al. (2008). 

The largest galaxy (G3) is chosen to have a stellar mass 
5 X 10^° Mq, and the smaller galaxies are chosen to have 
stellar masses 1.5 x 10^° M© (G2), 5 x 10^ M© (Gl), and 
1 x 10^ Mq (GO), spanning a factor of 50 in stellar mgiss. The 
stellar half-light radii are from the stellar mass-size relation 



of Shen et al. (2003). The bulge-to-disc ratios are taken from 
de Jong (1996) and used to determine the stellar disc and 
bulge masses and scalelengths. The gas fractions and masses 
are determined from the gas mass - stellar mass scaling re- 
lation from Bell et al. (2003). The gas scalelength is as- 
sumed to be three times the stellar disc scalelength (Broeils 
& Woerden 1994). We adopt NFW dark matter halo pro- 
files selected such that the rotation curves lie on the bary- 
onic TuUy-Fisher relation (Bell & de Jong 2001; Bell et al. 
2003). These models do not include adiabatic contraction. 
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Figure 3. u — r—z composite colour images including dust extinction for the 39:1 baryonic mass ratio tilted prograde-prograde GSGOPt 
simulation. The GO secondary galaxy is stripped and fades considerably, and its position is circled. GO is almost completely disrupted 
by the fourth pass. The images show the initial galaxies [1], the first pass [2], the maximal separation after the first pass [3], the second 
pass [4], the third pass [5], the fourth pass [6], and the last timestep calculated [7]. The viewing angles and image scales are the same as 
Figure [T] 



Table 2. Merger Mass Ratios 



Primary 


Satellite 


Total 


Stellar 


Baryonic 


G3 


G3 


1:1 


1:1 


1:1 


G3 


G2 


2.3:1 


3.3:1 


3.1:1 


G3 


Gl 


5.8:1 


10.0:1 


8.9:1 


G3 


GO 


22.7:1 


50.0:1 


38.9:1 


G2 


G2 


1:1 


1:1 


1:1 


G2 


Gl 


2.6:1 


3.0:1 


2.8:1 


G2 


GO 


10.0:1 


15.0:1 


12.4:1 


Gl 


Gl 


1:1 


1:1 


1:1 


Gl 


GO 


3.9:1 


5.0:1 


4.4:1 



The total mass-to-light ratio is assumed to vary with mass 
such that lower mass galaxy model have higher mass-to-light 
ratios. 

2.4 Merger Parameters 

The new galaxy merger simulations presented here are merg- 
ers of discs spanning a range of mass ratios and merger or- 
bits and orientations. Mergers of equal-mass gas-rich merg- 
ers were previously presented in Paper 1. The total, bary- 
onic, and stellar mass ratios are given in Table [J] We will 
generally refer the baryonic mass ratios of the merger, but 
note that the stellar /total mass ratios are higher/lower than 
the baryonic mass ratios. These ratios range from 1:1 — 23:1 
for total mass ratio, 1:1 — 39:1 for baryonic mass ratio, and 
1:1-50:1 for stellar mass ratio. The G3G3, G2G2, GlGl 



mergers have 1:1 mass ratios. The G3G2 and G2G1 mergers 
have baryonic mass ratios of ~ 3:1 and hence are more rep- 
resentative of the typical major merger than rare equal-mass 
mergers. The G3G1, G3G0, G2G0, and GIGO mergers have 
baryonic mass ratios of 9:1, 39:1, 12:1, and 4:1 respectively, 
and are considered to be minor mergers. 

Each permutation of primary-satellite merger has been 
simulated with a tilted prograde-prograde orientation with 
an initial sub-parabolic orbit (labeled Ft) with a pericen- 
tric radius Rpsri = 13.6 kpc. Additional simulations of 
major G3G2 mergers and minor G3G1 mergers have been 
performed with a wide range of orbital parameters (Table 
[3|. These include tilted prograde-prograde orientations on 
sub-parabolic orbit with two additional pericentric distances 
(P_: Rperi = 6.8 kpc, P+: Rperi = 27.2 kpc), pure prograde- 
prograde orientation (P), pure and tilted polar orientations 
(PI, Pit), primary retrograde — satellite prograde orientation 
(RP), primary prograde — satellite retrograde orientation 
(PR), and circular orbit (C; see Table|3}. All of simulations 
initially have slightly sub-parabolic orbits with an eccentric- 
ity e = 0.95, except for the circular orbits simulations, which 
by definition have e = 0. All of the parabolic orbit simula- 
tions start with an initial separation of 250 kpc; the circular 
orbit simulations start with an initial separation of 120 kpc. 

For comparison, we have also calculated the morpholo- 
gies for an isolated G3 galaxy. The isolated G3 simulation 
was run for 6 Gyr. Because no additional gas is provided, 
the isolated G3 galaxy consumes its initial gas reservoir and 
fades as the simulation progresses. 
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Table 3. Merger Simulation Parameters 



Simulation 










R d 




(deg) 


(deg) 


(deg) 


(deg) 




(kpc) 




Equal-mass prog 


rade-prograde merj 


;ers 




G3G3Pi^ 


-30 





30 


60 


0.95 


13.6 


G2G2Pj 


-30 





30 


60 


0.95 


Q Q 
O.O 


GlGlPf 


-30 





30 


60 


0.95 


3.0 




Unequal-mass pro 


grade-prog 


;rade mergers 




G3G2Pt 


-30 





30 


60 


0.95 


1 Q R 
io.O 


GSGlPt 


-30 





30 


60 


0.95 


iO.D 


G3G0Pt 


-30 





30 


60 


0.95 


io.O 


G2GlPt 


-30 





30 


60 


0.95 


O.oZ 


G2G0Pt 


-30 





30 


60 


0.95 


O.oZ 


GlGOPt 


-30 





30 


60 


0.95 


z.yo 




3:1 baryonic : 


mass ratio 


mergers 






G3G2P 








30 


60 


0.95 


13.6 


G3G2P_ 


-30 





30 


60 


0.95 


6.8 


G3G2P+ 


-30 





30 


60 


0.95 


27.2 


G3G2Plt 


-60 





30 


60 


0.95 


13.6 


G3G2P1 


-90 





30 


60 


0.95 


13.6 


G3G2RP 


180 





30 


60 


0.95 


13.6 


G3G2PR 


-30 





150 





0.95 


13.6 


G3G2C 


-30 





30 


60 


0. 




9:1 baryonic mass ratio mergers 


G3G1P 








30 


60 


0.95 


13.6 


G3G1P- 


-30 





30 


60 


0.95 


6.8 


G3G1P+ 


-30 





30 


60 


0.95 


27.2 


G3GlPlt 


-60 





30 


60 


0.95 


13.6 


G3G1P1 


-90 





30 


60 


0.95 


13.6 


G3G1RP 


180 





30 


60 


0.95 


13.6 


G3G1PR 


-30 





150 





0.95 


13.6 


G3G1C 


-30 





30 


60 


0. 





" Initial orientation of galaxy 1 with respect to the plane of the 

orbit in spherical coordinates, where 9 = arctan( v^^ ) and 
4) = arctan(^) 

^ Initial orientation of galaxy 2 

^ Eccentricity of the orbit, where a parabolic orbit has e = 1. 

Pericentric distance of the initial orbit 
^ In Paper 1, the equal-mass tilted prograde-prograde mergers 
are named G3PP, G2PP, and GIPP 



3 IMAGE ANALYSIS 

We replicate the observations and measurements of real 
galaxy mergers as closely as possible. We focus on rest-frame 
g morphologies for purposes of this paper, as these simula- 
tions can be used to calibrate the morphologies of galaxies 
currently observed from the ground and with the Hubble 
Space Telescope in optical and near-infrared wavelengths at 
< 2 < 3. In the following section we briefly describe how 
the simulated g images are degraded and analysed to match 
real galaxy morphology measurements; a more detailed dis- 
cussion may be found in Paperl. 



3.1 Image degradation 

The images are produced by SUNRISE for each simulation 
for 11 isotropically positioned viewpoints as a function of 
time from ~ 0.5 Gyr prior to the first pass to > 1 Gyr after 
the final coalescence in ~ 30 — 250 Myr timesteps depend- 
ing on the merger state, up to a maximum runtime of 6 
Gyr. Only the G3G1PR and G3G1C simulations have not 
coalesced by this time. The nuclei of the G3G0Pt simulation 
don't truely merge, but the GO satellite is almost completely 
disrupted by its fourth pass. The field of view of the output 
images ranges from 200 kpc during the initial stages and 
period of maximal separation to 100 kpc during the first 
pass, second pass, final merger and post-merger stages. The 
intrinsic resolution of the output SUNRISE p-band images 
is 333 pc per pixel. 

The images output by SUNRISE have no background 
sky noise and no seeing effects, although they do have parti- 
cle noise and Monte Carlo Poisson noise. We degrade these 
images to simulate real data, but do not attempt to mimic a 
particular set of galaxy survey observations. We re-bin the 
images to 105 pc per pixel and convolve the images with a 
Gaussian function with a FWHM = 400 pc. This was done 
to simulate the effect of seeing but maintain as high spatial 
resolution as possible. The values where chosen to match 
the Sloan Digital Sky Survey (Abazajian et al. 2003) with 
1.5" seeing, 0.396" per pixel plate scale for a galaxy at a dis- 
tance such that 1.5" ~ 400 pc. We also add random Poisson 
noise background to simulate sky noise but scale this noise 
to maintain a high signal-to-noise for the primary galaxies 
(> 20 per pixel within the Petrosian radius). For the G3G0 
simulation, the GO satellite galaxy fades substantially be- 
cause of tidal stripping and declining star- formation, and 
has average signal-to-noise per pixel < 3 after the third pass. 



3.2 Morphology Measurements 

Each image is run through an automated galaxy detection 
algorithm integrated into our IDL code. If the centres of the 
merging galaxies are less than 10 kpc apart, they are gen- 
erally detected as a single object. If 2 distinct galaxies are 
detected, the detection segmentation maps are used to mask 
out the other galaxy while each galaxys morphology is mea- 
sured. The projected separation Rproj is measured when two 
galaxies are detected. For each detected object, we calculate 
the Petrosian radii in circular and elliptical apertures, con- 
centration C , 180 degree rotational asymmetry A, the Gini 
coefficient G, and the second-order moment of the bright- 
est 20% of the light M20. (Please refer to Lotz, Primack, 
& Madau 2004, Conselice 2003, and Paper 1 for detailed 
definitions) . 

In Paper 1, we studied the effect of numerical resolution 
on the simulation morphologies. We found small but signif- 
icant differences in the average M20 and A values when we 
compared the standard numerical resolution simulation to a 
simulation with 10 times as many particles, and corrected for 
these offsets. Because the numerical resolution of the simu- 
lations presented here are similar to the standard resolution 
simulations in Paper 1 (~ 10^ particles), we apply the same 
correction of SM2Q — —0.157 and SA = —0.115 to the values 
in this work. 
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Figure 4. Time v. projected separation (Rpro]), logio(SFR), and morphology (A, G, M20, C) for G3G3Pt, G3G2Pt, GSGlPt, GSGOPt 
mergers and isolated G3 galaxy. The satellite galaxies are shown as crosses, and the primary galaxy/merger is shown as diamonds. 
The merger stages are colour coded with initial galaxies:red, the first pass:green, the maximal separationiblue, second pass (and third 
pass): purple, final mergeriorange, post-merger:magenta, and merger remnant at > 1 Gyr after coalescence of the nuclei as cyan. The 
isolated G3 points are colour-coded to match the G3G0Pt simulation. The simulations show peaks in star-formation and morphological 
disturbances at the first pass and final merger, with short-lived or no asymmetries for the minor mergers G3G1 and G3G0. 



3.3 Merger Classification and Time-Scales 

Lotz et al. (2004) found that local ultra-luminous galaxies 
visually classified as mergers could be distinguished from the 
sequence of normal Hubble type galaxies with 

G > -0.115 A/ao +0.384 (1) 

or 

G > -0.4 yl -I- 0.66 or A > 0.4 (2) 

Asymmetry alone is also often used to classify merger 
candidates. The calibration of local mergers by Conselice 
(2003) finds the following merger criterion: 

A > 0.35 (3) 



Galaxies at higher redshift cannot be imaged at as high 
spatial resolution as local galaxies even when observed with 
HST. The measured morphologies of galaxies at 2; > 0.25 
imaged with HST will have non-negligible biases as a result 
of this lower spatial resolution (Lotz et al. 04). See Paper 1 
for a discussion of how the morphologies and merger criteria 
change for z ^ 1 HST resolution, and how the time-scale 
computed here may be applied to HST data. 

Close kinematic pairs are also probable pre-merger sys- 
tems. We assume h — 0.7 and we estimate the time-scales 
during which merging galaxies can be found as separate ob- 
jects within 5 < Rproj < 20, 10 < Rproj < 30, 10 < Rproj < 
50, and 10 < Rproj < 100 kpc. The simulated merging 
galaxies always have relative velocities < 500 km s~^. 

We calculate each simulation's average observability 
time-scale for the G — M20, G — A, and A criteria given 



© 2009 RAS, MNRAS 000, [T]-?? 



8 J.M. Lotz et al 




fim Tffijr) TtV) MOy^} TfOyi-) 



Figure 5. Time v. projected separation (Rproj), logio{SFR), and morpliology (A, G, A/20, C") for lower-mass G2G2Pt, G2GlPt, G2G0Pt, 
GlGlPt and GlGOPt mergers. The satellite galaxies are shown as crosses, and the primary galaxy/merger is shown as diamonds. The 
merger stages are colour coded as in the previous figure. As for the GSGlPt minor merger, the minor mergers G2G0Pt and GlGOPt 
appear asymmetric for a much shorter period of time than the major mergers G2G2Pt, G2GlPt, and GlGlPj. 



above by averaging the results of the 11 isotropic view- 
points. Because we wish to determine the number density of 
merger events rather than the number of galaxies undergoing 
a merger, galaxies that have not yet merged but identified 
morphologically as merger candidates are weighted accord- 
ingly. The time that each pre-merger galaxy is morphologi- 
cally disturbed is added (not averaged) to the time that the 
post-merger system appears disturbed. No such weighting is 
done for the close pair time-scales as this factor is generally 
included in the merger rate calculation (e.g. Patton et al. 
2000). 



4 RESULTS 

4.1 Evolution with Merger Stage 

The time evolution of each simulation depends on the mass 
ratios, masses, and orbits of the galaxies. Therefore it is use- 
ful to compare the behavior of each simulation as a function 
of merger stage rather than time. We define seven different 
merger stages based on the positions of the galaxy nuclei 
in physical space, which are colour-coded in Figures U}- [S] 
These are the initial encounter (red), the first pass (green), 
maximal separation immediately after the first pass (blue), 
the second pass or final approach (purple), the final merger 
(orange), the post-merger observed 0.5—1 Gyr after the final 
merger (magenta), and the remnant observed > 1 Gyr after 
the final merger (cyan). 

In Figures |4] and [5] we examine the projected separa- 
tions Rproj, measured morphologies (G, Af20, C, and A), 
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Figure 6. G-hho, G-A, and C-Afor for G3G3Pi, G3G2Pt, GSGlPt, G3G0Pt mergers and isolated G3 galaxy. The satellite galaxies 
are shown as crosses, and the primary galaxy/merger is shown as diamonds. The merger stages are colour-coded as in the previous plots. 
The isolated G3 galaxy is colour-coded to match the stages of the G3G0 simulation. The simulations show peaks in morphological 
disturbances at the first pass (green points), second pass (purple points), and final merger (orange points), with short-lived or no 
asymmetries for the minor mergers G3GlPt and G3G0Pt . 



and star-formation rate per object as a function of time and 
merger stage for the tilted prograde-prograde simulations, 
spanning a range of baryonic mass ratio of 1:1 — 39:1 for 
the G3 simulations (Figure IH and 1:1 — 12:1 for the G2 
and Gl simulations (Figure [SJ. The scatter in Rproj and 
morphologies at each timestep is the variation in the merger 
appearance with the 11 viewing angles. The minor mergers 
GSGlPt and G2G0Pt take significantly longer to complete 
than the major mergers, and experience a third and final 
pass before the galaxies coalesce. In the GSGOPt, the GO 
satellite is disrupted and becomes too faint to be detected 
above the noise in our simulated images after the third pass 
so the projected separation falls to zero (Figure [3]). 

The simulations generally show similar correlations of 
star-formation rate and morphology with merger-stage. The 
initial segmentation maps computed to identify each galaxy 
are used to compute the total star-formation rate for each 
galaxy at each timestep and camera. The major mergers 
experience a burst of star-formation at the first pass and 
a few hundred Myr after the galaxies merge as their nuclei 
coalesce. As we found in Paperl, this peak in star-formation 
rate at the final merger occurs a few hundred Myr after the 
strongest morphological disturbances in A, C, and M20. 

The minor mergers also have enhanced star-formation 
rates at the close passages. However, the minor mergers' 
burst efficiencies, or stars formed in excess of what is ex- 
pected for isolated galaxies, is lower than for major merg- 



ers (Cox et al. 2008; see also Johansson et al. 2009). The 
G3GlPt and G2G0Pt mergers do not show strong starbursts 
at the final merger. This is because the systems have already 
consumed most of their gas by the time of the final merger 
~ 2 — 4 Gyr after the first encounter. The G3G0Pt simula- 
tion shows the strongest decline in the star-formation rate of 
all the G3 merger simulations. The G3G0 has the least gas 
available for star-formation because of the low gas mass of 
the GO satellite, and the evolution of its star-formation rate 
with time is very similar to the isolated G3 galaxy (Figure 

Except for the G3G0Pt simulation, all of the mergers 
show peaks in A and M20 and minima in C during the close 
passages and just before the final merger, as was also found 
for the equal-mass mergers in Paper 1. However, mass ra- 
tio affects the strength of these morphological disturbances. 
In Figure |4l we find that the major mergers G3G3Pt and 
G3G2Pt have roughly similar offsets in A and C at the 
first pass and final approach/merger. The minor merger 
G3GlPt has much weaker disturbances in A and C at the 
first pass and final approach/merger, while G3G0Pt shows 
only a small increase in A at the first pass. The increase in 
AI20, on the other hand, is equally strong for the G3G3Pt, 
G3G2Pt, and G3GlPt mergers. In Figure O we see that the 
lower mass G2 and Gl mergers show similar dependence 
on mass ratio. The major mergers G2G2Pt, G2GlPi, and 
GlGlPt have similarly strong peaks in A and minima in C 
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Figure 7. G- M20, G - A, and C-Aior lower-mass G2G2Pt, G2GlPt, G2G0Pt, GlGlPt and GlGOPt mergers. The satellite galaxies 
are shown as crosses, and the primary galaxy/merger is shown as diamonds. The merger stages are colour coded as in the previous figure. 
As for the GSGlPt minor merger, the minor mergers G2G0Pt and GlGOPt are asymmetric for a much shorter period of time than the 
major mergers G2G2Pt, G2GlPt, and GlGlPt. 
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Figure 8. Time-scales for morphological disturbances in G — M20, G — A, and A v. baryonic mass ratio. The black crosses show the 
time-scales for each viewing angle, and the red diamonds are the viewing-angle averaged time-scales. The time-scales for observing high 
asymmetry is a strong fonction of mass ratio, while the time-scales for observing objects with high G — M20 values appears to be 
independent of mass ratio below baryonic mass ratio ~ 10. Therefore G — M20, unlike A, will identify minor mergers with baryonic mass 
ratios < 9:1 as well as major mergers. 



at the first pass and final approach/merger, while the minor 
mergers G2G0Pt and GlGOPt have weaker disturbances in 
A and C. 

Dust extinction is a function of the metallicity and gas- 
fraction of the merging system, and has the strongest effect 
on G (Paper 1). For the more massive G3 and G2 mergers 



with higher gas metallicity, G shows significant scatter with 
viewing angle between the first pass and final merger as a 
result of dust extinction along certain lines of sight. The 
lowest mass Gl mergers show higher G values between the 
first pass and final merger because this simulation has the 
lowest metallicity and dust content. 
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Table 4. Morphological Disturbance Time-Scales 



Simulation T(G - M20) T(G - A) T{A) 
(Gyr) (Gyr) (Gyr) 



Equal-mass prograde-prograde mergers 



G3G3Pt 


0.16± 0.07 


0.31± 0.08 


0.23± 0.11 


G2G2Pt 


0.22± 0.14 


0.31± 0.15 


0.25± 0.19 


GlGlPt 


0.23± 0.04 


0.35± 0.14 


0.30± 0.16 


Unequal-mass prograde-prograde mer 


■gers 


G3G2Pt 


0.25± 0.08 


0.30± 0.13 


0.24± 0.11 


GSGlPt 


0.36± 0.15 


0.27± 0.13 


0.03± 0.03 


GSGOPt 


0.03± 0.03 


0.02± 0.03 


O.OOit 0.00 


G2GlPt 


0.17± 0.06 


0.30± 0.09 


0.26± 0.10 


G2G0Pt 


O.lOit 0.08 


0.06± 0.05 


O.Olib 0.02 


GlGOPt 


O.llit 0.06 


0.22± 0.10 


0.14± 0.07 


3:1 baryonic mass ratio mergers 


G3G2P 


0.20± 0.10 


0.28± 0.06 


0.25± 0.06 


G3G2P_ 


0.32± 0.10 


0.39± 0.12 


0.30± 0.12 


G3G2P+ 


0.20± 0.14 


0.27± 0.10 


0.18± 0.11 


G3G2Plt 


0.30± 0.14 


0.36± 0.09 


0.25± 0.09 


G3G2P1 


0.34± 0.11 


0.47± 0.14 


0.37± 0.10 


G3G2RP 


0.38± 0.14 


0.43± 0.11 


0.35± 0.11 


G3G2PR 


0.28± 0.11 


0.44± 0.17 


0.29± 0.07 


G3G2C 


0.31± 0.17 


0.40± 0.20 


0.36± 0.17 




9:1 baryonic 1 


mass ratio mergers 




G3G1P 


0.24± 0.11 


0.21± 0.13 


0.03± 0.03 


G3G1P- 


0.26± 0.07 


0.26± 0.14 


0.06± 0.04 


G3G1P+ 


0.46± 0.26 


0.24± 0.06 


0.02± 0.03 


G3GlPlt 


0.37± 0.13 


0.26± 0.15 


O.Olit 0.04 


G3G1P1 


0.32± 0.08 


0.28± 0.14 


0.06± 0.05 


G3G1RP 


0.22± 0.07 


O.llit 0.10 


0.06± 0.05 


GSGIPR^ 


> 0.26± 0.13 


> 0.08± 0.08 


> 0.03± 0.07 


GSGIC" 


> 0.05± 0.06 


> O.OOit 0.01 


> O.OOit 0.00 



" The simulation does not merge within 6 Gyr. 

As in Paper 1, we find that the major merger remnants 
observed > 1 Gyr after the merger have G, M20, C, A, and 
star-formation rates more similar to local early-type spiral 

galaxies (G ~ 0.55, A/20 2.0, C ~ 3.5, A ~ 0.05) than 

local elliptical/so galaxies (G ~ 0.60, M20 ~ -2.5, G - 4.5, 
A ~ 0.0). These merger remnants retain significant amounts 
of gas and reform small-mass star- forming discs. Despite 
their massive bulges, these low-mass disc components are 
bright enough in p-band to affect the remnants' quantita- 
tive morphologies (see also Barnes 2002, Springel et al. 2005, 
Naab et al. 2006, Robertson et al. 2006). We note that our 
simulations do not include any form of AGN feedback, which 
could remove this remaining gas and prevent the formation 
of discs large enough to affect the remnant morphologies 
(e.g. di Matteo et al. 2005), although AGN feedback may 
not be sufficient to quench star-formation in unequal-mass 
mergers (Johansson et al. 2009). The supernovae feedback 
efficiency and merger geometries can also effect the amount 
and spatial distribution of gas in the merger remnant (e.g. 
Springel et al. 2005, Cox et al. 2006, Paper 1). 

The minor merger remnants, on the other hand, do have 
quantitative morphologies similar to ellipticals/SO galaxies. 



But these simulations have been run for several Gyr longer 
than the major merger simulations without any additional 
gas supplied. Visual inspection of the images reveals that 
the minor merger remnants have compact nuclei and red 
discs rather than massive bulge components (Figures (2] |3]) . 
Thus the difference in the morphologies for major and mi- 
nor merger remnants is likely the result of longer simulation 
runtimes/greater gas consumption for the minor merger sim- 
ulations. We find that an isolated G3 galaxy simulation run 
for 6 Gyr with no additional gas supply results in an ob- 
ject very similar in appearance and quantitative morphol- 
ogy to the minor merger remnants (Figure I?}. Therefore we 
conclude that minor gas-rich mergers have a minimal effect 
on the morphologies of their merger remnants. The optical 
properties of the merger remnants will be explored in more 
detail in a future paper. 

4.2 Morphology Time-Scales 

We calculate the viewing-angle averaged time each merger 
simulation exhibits morphologies which meet the merger cri- 
teria for G - M20, G-A, and A (Table i]). We find that the 
mass ratio is an important factor in determining how long or 
whether a merger may be identified as a merger in G — M20, 
G — A, OT A, but the merging galaxy orientations and orbits 
are not. 

The merger simulations clearly show different behav- 
iors in G — A/20, G — A, and yl as a function of mass 
ratio. We plot G - A/20, G-A, and G - A for the 
GSGSPt, G3G2Pt, GSGIP*, GSGOPt mergers (Figure [5} 
and G2G2Pt, G2GlPt, G2G0Pt, GlGlPt, and GlGOPt 
mergers (Figure (Tj . We find that a significant number of 
timesteps lie in the merger region of the G — A/20 plots 
for all the merger simulations except for the highest mass 
ratio GSGOPt simulation. On the other hand, only the 
major mergers (G3G3Pt, G3G2Pi, G2G2Pt, G2GlPt, and 
GlGlPt) have a large number of timesteps in the merger 
regions of the G — A and A plots. The isolated G3 simula- 
tion never shows morphologies that meet the merger criteria 
(right panels, Figure 

In Figure [S] we plot the observability time-scales v. 
baryonic mass ratio for the tilted prograde-prograde merg- 
ers. We find that G — A/20 shows no significant decrease in its 
observability time-scale for baryonic mass ratios < 9:1, and 
is as likely to find 9:1 minor mergers as 1:1 major mergers. 
The asymmetry time-scales are stronger function of mass 
ratio, and A > 0.35 does not detect mergers with baryonic 
mass ratios > 5:1. 

We visually inspected the merger images to better un- 
derstand why G — A/20 and A show these different be- 
haviours. In Figure [S] we show the g-hand images of tilted 
prograde-prograde mergers viewed just before the final 
merger (or at the second pass for G3G0Pt). A/20 is a mea- 
sure of the second-order moment of the brightest 20% of an 
object's fiux, normalized by the total second-order moment. 
The brightest 20% of the light is shown in Figure [9] by the 
inner red contours, and the outer black contour shows the 
segmentation map used to compute G and A/20 • The nucleus 
of the merging satellite is among the brightest 20% of the 
light at this merger stage for all the mergers except for the 
highest mass ratio merger G3G0Pt. Therefore, G — A/20 will 
detect mergers when the nucleus of the merging satellite is 
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among the brightest 20% of the system's pixels. Asymmetry, 
on the other hand, is a measure of the rotational asymmetry 
of the system, and depends on the sum of the absolute differ- 
ence between the original image and the image rotated 180 
degrees about the center of minimum asymmetry. Therefore 
A depends more strongly on the relative brightnesses of the 
nuclei, and is largest when the nuclei are of similar luminos- 
ity. 

As we found for equal-mass mergers in Paper 1, the 
orbits and relative orientations do not have a significant ef- 
fect on the time-scales for identifying unequal-mass galaxy 
mergers for any of the morphology criteria. In Figure [TOl we 
plot the observability time-scales for G — M20, G — A, and A 
for the G3G2 and G3G1 mergers with nine different orbits 
and orientations. The G3G1 circular orbit (C) and G3 pro- 
grade — Gl retrograde parabolic orbit (PR) simulations do 
not merge within simulation runtimes of 6 Gyr, and hence 
the plotted time-scales are upper limits (arrows) . The mean 
observability time-scales for each orbit/orientation (red di- 
amonds) have less variation than the scatter due to viewing 
angle for a given orbit/orientation (black crosses). 

4.3 Close Pair Separation Time-Scales 

The time-scales for identifying the simulations as close pairs 
with projected separations 5 < Rproj < 20, 10 < Rproj < 
30, 10 < Rproj < 50, and 10 < Rproj < 100 kpc h''^ in 
Tableland Figure [TT] As expected from dynamical friction 
arguments, the close pair time-scales are longer for higher 
mass ratio mergers. The relative orbits and orientations are 
less important for the observability time-scales (Figure [TT|) . 
except when the galaxies are at large separations on circular 
orbits or start with larger pericentric distances. 

We compare our results to the merger time-scales of 
Boylan-Kolchin, Ma & Quataert (2008) for GADGET-2 
mergers of dark matter haloes, and of Kitzbichler & White 
(2008) derived from the Millenium Simulation semi-analytic 
model. The dynamical decay time-scale for the G3-G2 
mergers are ~ 1-2 Gyr shorter than Boylan-Kolchin et al. 
(2008) dynamical friction time-scale for the dark matter halo 
merger of similar total mass, slightly higher total mass ratio 
(3:1), initial separation, and similar orbits, possibly because 
of the effect of baryons. The minor G3-G1 merger time-scales 
for sub-parabolic orbits (~ 4 Gyr) appear to be more con- 
sistent with Boylan-Kolchin et al. (2008) dynamical friction 
time-scale for mergers with similar orbits and mass ratios 
(~ 3.5-4.5 Gyr). 

The time-scale for a G3-G2 merger on a circular or- 
bit to appear as a close pair with projected separation 
10 < Rproj < 50 kpc (1.2 ± 0.6 Gyr) is consistent with 
Kitzbichler & White (2008) calculations of the close pair 
time-scale for a 5 x 10^" M0 stellar mass galaxy pair ob- 
served with Rproj < 50 kpc (1.4 Gyr). The Kitzbichler 
& White (2008) close pair time-scales are based on the Cro- 
ton et al. (2006) semi-analytic model for galaxies in the dark 
matter haloes of the Millennium Simulation (Springel et al. 
2005), and assume the Binney & Tremaine (1987) dynam- 
ical friction time-scale for objects on a circular orbit when 
the subhaloes hosting galaxies can no longer be resolved in 
the Millennium Simulation. When compared to our circular 
orbit simulations, the majority of parabolic orbits simula- 
tions yield ~ 15-30% shorter close pair time-scales. At larger 
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Figure 9. The brightest 20% of light is marked in red contours 
for G3G3Pt, G3G2Pt, GSGlPt, and GSGOPt mergers just before 
the final coalesence (except for GSGOPt, which is shown at the 
second pass and does not merge). The sateUite nuclei are among 
the brightest regions for mergers with baryonic mass ratio < 9:1 
(G3G1), but are too faint to be detected by G— A/20 at baryonic 
mass ratio 39:1 (i.e. G3G0). 

projected separations {Rproj < 100 kpc h~^), we find the 
parabolic orbit simulations have ~ 40% shorter time-scales 
than the circular orbit simulations. This implies that the 
assumption of a circularized orbit may result in ~ 15-30% 
over-estimation of the typical similar mass close pair time- 
scale at Rproj < 50 kpc and ~ 40% over-estimation of 
typical close pair time-scales at Rproj < 100 kpc . 

We have compared the time-scales for the true sepa- 
rations of the mergers to those derived for the projected 
simulations. In Figure 1121 we show that the time-scales for 
observing a galaxy pair at a given projected separations are 
~ twice as long as the time-scale for the galaxy pair at the 
same separation in physical space. This is because projec- 
tion efi'ects will always make galaxy pairs appear closer than 
they truly are. 



5 SUMMARY AND IMPLICATIONS 

We have analysed the g-band quantitative morphologies 
and projected pair separations for a large suite of GAD- 
GET/SUNRISE simulations of unequal-mass disc galaxy 
mergers. These merger simulations span a range in baryonic 
mass ratio from 1:1 to 39:1, and include a range of orienta- 
tions and orbital parameters for the merging galaxies. The 
initial galaxies have baryonic gas fractions matched to local 
disc galaxies. We determine the observability time-scales for 
identifying these simulated mergers using the quantitative 
morphology classification in G — M20, G — A, and A and 
as close pairs with projected separations 5 < Rproj < 20, 
10 < Rproj < 30, 10 < Rproj < 50, and 10 < Rproj < 100 
kpc . Our main conclusions are as follows: 
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Figure 10. Time-scales for morphological disturbances in G — A/20, G — A, and A v. orientation and orbits for the 3:1 baryonic mass 
ratio G3G2 and 9:1 baryonic mass ratio G3G1 mergers (see Table The time-scales do not depend strongly on the merger orientations 
or orbits. The black crosses show the time-scale for each viewing angle, and the red diamonds are the viewing-angle averaged time-scales. 
The G3G1PR and G3G1C simulations do not merge within 6 Gyr, hence the time-scales shown are lower limits (arrows). 
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Figure 11. Time-scales for projected separations of 5 < Rproj < 20, 10 < Rproj < 30, 10 < Rproj < 50, and 10 < Rproj < 100 
kpc V. orientation/orbits for baryonic mass ratio 3:1 G3G2 mergers (top) and baryonic mass ratio 9:1 G3G1 mergers (bottom). The 
time-scales are significantly longer for circular orbit mergers observed at large projected separations or for high mass ratios. Symbols 
and orbits are the same as previous figure. 
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Figure 12. Left: The close pair time-scales calculated for real separations v. real maximum separation for all G3-G2 merger simulations. 
Center: The average close pair time-scales calculated for projected separation v. projected maximum separation for all G3-G2 merger 
simulations. The error-bars are the standard deviations with viewing angle. Right: The average projected separation time-scale divided by 
real separation time-scale v. maximum separation. The projected separation time-scales average about twice that of the decay time-scales 
at that separation in real space because projection effects make galaxies appear closer than they truly are. 



Table 5. Close Pair Time-Scales 



Simulation 


T(5 < Rproj < 20) T(10 < Rproi < 30) 


T{10 < Rpro3 < 50) 


T (10 < Rpro] < 100) 




(Gyr) 


(Gyr) 


(Gyr) 


(Gyr) 






Equal-mass prograde-prog 


radc mergers 




G3G3Pt 


0.39± 0.30 


0.72± 0.39 


1.21± 0.38 


1.85± 0.13 


G2G2Pi^ 


0.43± 0.19 


0.60± 0.16 


0.72± 0.20 




GlGlPt" 


0.59± 0.12 


0.52± 0.16 


0.64± 0.19 








Unequal-mass prograde-pro 


grade mergers 




G3G2Pt 


0.26± 0.15 


0.59± 0.44 


1.15± 0.36 


1.98± 0.16 


G3GlPt 


0.44± 0.28 


0.96± 0.47 


1.67± 0.62 


2.87± 0.31 


G3G0Pt 


0.58± 0.42 


1.23± 0.53 


2.75± 0.54 


3.44± 0.34 


G2GlPt" 


0.48± 0.18 


0.67± 0.18 


0.76± 0.22 




G2G0Pi" 


l.Olit 0.32 


1.17± 0.35 


1.64± 0.29 




GlGOPt" 


0.79± 0.32 


0.78± 0.31 






3:1 baryonic mass ratio mergers 


G3G2P 


0.34± 0.29 


0.62± 0.48 


1.24± 0.42 


1.92± 0.18 


G3G2P_ 


0.28± 0.15 


0.44± 0.16 


0.91± 0.23 


1.36± 0.28 


G3G2P+ 


0.30± 0.17 


O.SOit 0.51 


1.41± 0.74 


2.96± 0.33 


G3G2Plt 


0.27± 0.17 


0.62± 0.47 


1.18± 0.44 


1.97± 0.14 


G3G2P1 


0.29± 0.20 


0.58± 0.42 


1.15± 0.40 


2.00± 0.15 


G3G2RP 


0.30± 0.32 


0.59± 0.48 


1.19± 0.46 


1.96± 0.11 


G3G2PR 


0.24± 0.27 


0.57± 0.48 


1.21± 0.51 


1.92± 0.13 


G3G2C 


0.32± 1.38 


0.82± 0.31 


1.42± 0.64 


3.20± 0.84 


9:1 baryonic mass ratio mergers 


G3G1P 


0.49± 0.24 


0.96± 0.49 


1.58± 0.62 


2.80± 0.34 


G3G1P_ 


0.33± 0.25 


0.62± 0.26 


1.25± 0.36 


1.80± 0.33 


G3G1P-I- 


0.57± 0.49 


1.12± 0.69 


2.22± 0.98 


4.00± 1.14 


G3GlPlt 


0.44± 0.28 


l.Olit 0.49 


1.70± 0.74 


3.05± 0.36 


G3G1P1 


0.42± 0.26 


0.91± 0.48 


1.67± 0.73 


3.06± 0.34 


G3G1RP 


0.57± 0.39 


1.05± 0.44 


1.68± 0.61 


3.03± 0.37 


G3G1PR^ 


> 0.43± 0.25 


> l.Olit 0.41 


> 1.64± 0.63 


> 3.04± 0.27 


G3G1C''~ 


> 0.18± 0.31 


> 1.23± 1.16 


> 2.64± 1.16 


> 5.06± 0.95 


Time-scale 


is not calculated when the initial simulation separation is less than 


maximum projected se 



^ The simulation does not merge within 6 Gyr. 
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• Different morphological approaches identify mergers 
with different mass ratios. G — M20 is sensitive to mergers 
with baryonic mass ratios between f :f and at least 9:f , and 
does not show a correlation of observability time-scale with 
baryonic mass ratio less than 9:1. Typical G — M20 time- 
scales are 0.2 — 0.4 Gyr. Asymmetry, on the other hand, is 
much less likely to find minor mergers with baryonic mass 
ratios greater than 4:1 when gas- fractions typical of local 
galaxies are assumed. The asymmetry time-scales for major 
mergers of moderate gas- fraction discs are also 0.2 — 0.4 Gyr, 
and less than 0.06 Gyr for minor mergers of moderate gas- 
fraction. Hence, the G — M20 method detects both major and 
minor mergers, while the majority of asymmetric objects in 
the local universe will be major mergers. From visual inspec- 
tion of the simulation images, we find that G — A/20 detects 
mergers when double nuclei are enclosed in a common enve- 
lope. If the nucleus of the merging satellite galaxy is among 
the brightest 20% of the pixels, it will produce high M20 val- 
ues at close passages and final merger. A is more sensitive to 
the relative brightnesses of the nuclei, and is strongest when 
the merging nuclei are of similar luminosity. 

• The orbits and orientations have little effect on the 
overall time-scale for strong morphological disturbances, but 
do effect the time-scale for identifing galaxy pairs. The time- 
scales for close pairs with projected separations Rproj < 
50 kpc are ~ 15—30% longer for circular orbits than 
parabolic orbits. For close pairs with projected separations 
> 100 kpc , we find that circular orbits result in close 
pair time-scales ~ 40% greater than those found for mergers 
initially on parabolic orbits. 

• Our major merger remnants observed > 1 Gyr af- 
ter the final merger have quantitative morphologies and 
star-formation rates more like early-type spirals than ellip- 
tical/SO galaxies. We note that only supernova feedback is 
including in our simulations, and other quenching mecha- 
nisms (e.g. AGN feedback) may prevent the reformation of 
a disc and additional star-formation. We find that minor 
mergers have a minimal effect on the G, M20, C, and A 
values of their remnants. 

In Paper 1, we found that the various methods for find- 
ing galaxy mergers identify a given equal-mass merger at dif- 
ferent stages in the merger process and for varying periods of 
time. However, these different time-scales for major mergers 
are not sufficient to explain the observed differences between 
the fractions of galaxies in close pairs and the fraction of 
galaxies with visually or quantitatively disturbed morpholo- 
gies (e.g. Lotz et al. 2008a, Jogee et al. 2009). With this 
work, we now find that the two main quantitative morphol- 
ogy approaches - G — M20 and A - axe sensitive to different 
ranges in merger mass ratios. 

The range of merger mass ratio detected by the dif- 
ferent techniques has important implications for computing 
the galaxy merger rate. In a given galaxy survey, the galaxy 
merger candidates identified by G — M20 contain both minor 
mergers and major majors, while A is more likely to find 
major mergers (assuming local disc galaxy gas fractions). 
Because G — M20 identifies the merger stages when galax- 
ies exhibit double nuclei, these minor mergers may also be 
found by visual classification. Close kinematic pair studies 
typically limit the mass or fiux ratios of the paired galaxies 
to < 2:1 — 4:1, and hence select major merger events. Thus 
our results here may resolve many of the discrepancies in 



the merger fractions and merger rates found using different 
methods. 

Finally, we note that we have not explored the proper- 
ties of mixed-morphology mergers (e.g. disc-spheroid: Jo- 
hansson et al. 2009) or spheriod-spheroid mergers (e.g. 
Naab, Khochfar, & Burkert 2006; BeU et al. 2006; Boylan- 
Kolchin, Ma, & Quataert 2005), which are likely to give 
shorter morphological disturbance time-scales. We expect 
the time-scales derived in this paper to be upper limits to 
the morphological disturbance time-scales of gas-poor merg- 
ers. We explore the effect of gas fraction on the morphologies 
of merging discs in a companion paper to this work (Paper 
3), and find that like mass ratio, gas-fraction is an important 
parameter in the morphological disturbance time-scales and 
the calculation of the galaxy merger rate. 

This series of papers are the needed first steps for under- 
standing which parameters are important for the time-scales 
of morphological disturbances, how these relate to close pair 
time-scales, and for deriving a merger rate from the obser- 
vations of morphologically disturbed galaxies. Although the 
parameter space explored in these works is not exhaustive, 
we now understand which parameters are most important 
for the duration of morphological disturbances — mass ra- 
tio and gas-fraction. We will apply these merger observabil- 
ity time-scales to observations of the merger fraction using 
G — M20, A, and close pairs in a future paper. By comparing 
the minor -I- major merger frequency detected by G — M20 
to the major merger frequency detected by either A or close 
pairs, we will be able to place some of the first quantitative 
constraints on the minor merger rate. 
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